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ABSTRACT 


A proposed navigation ccmcept requires that a user measure the time-delay that 
satellite- emitted signals experience in traversing the distance between satellite 
and user. Simultaneous measurement of the propagation time from four differ- 
ent satellites permits the user to determine his position and clock bias if satel- 
lite ephemerides and signal propagation velocity are known. A pulse propagating 
throu^i the ionosphere is slowed down somewhat, giving an apparent range that 
is larger than the equivalent free space range. The difference between the ap- 
parent range and the true range, or the free space velocity and the true velocity, 
is the quantity of interest. Iliis quantity is directly proportional to the total 
electron content along the path of the propagating signal. Thus, if the total elec- 
tron ccmtent is known, or is measured, a perfect correction to ranging could be 
performed. 

Faraday polarization measurements are continuously being taken at Fort Mon- 
mouth, N. J. , using beacon emissions of the ATS-3 (137.35MHz) satellite, 
which is in a geostationary orbit. The polarization data yield electron content 
values up to ~ 1500 km since the measurement technique is based on the Faraday 
effect which is weighted ly the geomagnetic field. Day-to-day variability of the 
diurnal variation of total electroi; content values is present with differences of 
up to 50% or more not being unconimon. In addition, superposed on the overall 
diurnal variation are smaller scale variations of ~ 5 to 10% of the total content 
which are attributed to ionospheric density irregularities. 

Future experiments planned for the emissions of the forthcoming ATS-F will 
permit Faraday rotation, dispersive phase, and dispersive group delay measure- 
ments. The latter two will give the integrated electron density to geostationary 
altitudes while the first will give the density integrated to ~ 1500 km. The dif- 
ference-referred to as the exospheric content— will yield the currently unknown 
propagation time delay from ground to geostationary altitudes. 


INTRODUCTION 

A space-based radio navigation system could provide military and civilian users 
with precise three-dimensional position and velocity data. 


427 



The navigation signals cmitain data such as satellite identity, real-time ephem- 
erides information, system time and other data. To determine his position and 
velocity, the user cross-correlates the coded time signal received from the 
satellite with the same coded time signal generated in his receiver. The relative 
phase, or equivalently the time displacement between the user's receiver and 
the incoming code, determines the range to the satellite. Simultaneous meas- 
urment of such relative phases from four different satellites permits the user 
to determine his range and his clock bias with respect to the satellite's position 
and clock respectively. In addition to such range measurements, corresponding 
range-rates are measured by the carrier frequency Doppler shift of each signal 
and the corresponding motion of the particular satellite as described by the 
ephemerides data. 

The range from user to the various satellites is, of course, not the correct geo- 
metric range. A propagating navigation signal is slowed down the ionosphere 
1^ an amount proportional to the total electron content along its path. The elec- 
tron content may be measured in real time, provided the user has dual-frequency 
capabilities. However, substantial redtK:tion in the cost of user equipment could 
be realized if the navigation system uses only one frequency. In such a case, the 
ionospheric time delay will have to be determined through empirical modelling 
techniques, based on existing and future global electrcm content data, and will 
be transmitted to the user for correction via the navigating signal. To calculate 
the true range, one must determine the group velocity of the signal along the 
path. 


The transit time of a transionospheric signal from a satellite to a ground observer 
is: 


t = 



( 1 ) 


where ds is an element of distance along the signal's path and Vg and ng are the 
gi'oiQ) velocity of the signal and the group refraction index respectively, c is the 
speed of light, and O and S are the observer and satellite position respectively. 
In the high-frequency approximation, the group index of refraction is: 


40.3 N 

n_ = I + , (2) 
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where N is the electron density per meter cubed, and f is the operating frequency 
in Hz. Equation 1 becomes: 

1 

c 

The first term in Equation 3 represents the free space signal transit time, while 
the second term represents the signal’s excess time delay in the ionosphere. 

The excess time delay is inversely proportional to the frequency squared so ^t 
in the L band (the proposed navigation frequency is in the 1600 MHz band) it is 
very small. However, the precision req\dred by the system is such that the 
ionospheric time-delay has to be accoimted for. 

Thevariatlonof the time delay \vith electron content for 1.6 GHz, 400 MHz, and 
150 MHz is shown in Figure 1. At 1.6 GHz, the excess delay times are 0.524 nano- 
seconds and 52. 4 nanoseconds for a total content of 10 ‘^el/m^ and lO^^el/m^ respec- 
tively, which are the lower and upper bounds of the normal vertical ionospheric 
electron content. For an oblique signal path, the total content increases by an 
amount proportional to the additional path length (assuming spherical stratifica- 
tion of the ionospheric density distribution). At low elevations, i.e. below 10°, 
the slant total content exceeds the vertical content (i.e. elevation angle = 90°) 
by a factor of 3; while at 40° elevation, the slant content exceeds the vertical 
content by a factor of 1. 5. Thus, for the vertical content at the upper bound and 
a low-elevation signal path, the total time delay introduced by the ionosphere is 
--150 nanoseconds at 1.6 GHz. 


40.3 

ds+ I N 
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ds. 


(3) 


THE FARADAY ROTATION 


In the high-frequency and quasi-longitudinal approximations, the two magneto- 
ionic modes are nearly circularly polarized in opposite senses, thus a plane 
polarized wave traversing the ionosphere may be regarded as the vector sum of 
the ordinary and extraordinary components. Since the two components travel 
at different phase velocities, the plane of polarization rotates continually along 
the signal's path. The total rotation from the signal source to the observer is 
related to the total electron content by the expression: 



r ^ 

I f2 

•'n 



(B cos 6 secx) N dh, 


(4) 


where k = 2.36 x 10'^ , B is the local magnetic field flux density in gammas, 6 
is the angle between the radio wave normal and the magnetic field direction, and 
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Figure 1. Ionospheric Excess Time Delay vs. Integrated Electron Ck>ntent 
for 150 MHz, 400 MHz, and 1600 MHz Signals 



X is the angle between the wave normal and the vertical. Since B decreases 
inversely with the cube of the geocentric distance, and since the electron den- 
sity decreases exponentially with altitude above F 2 max (^?00km), the integral 
is heavily weighted near the earth and is considered to provide electron content 
values below ~1500km. 

The term M = B cos 6 sec x In Equation 4 may be taken out of the integral sign 
and replaced its value at a ’’mean" ionospheric altitude. Equation 4 then 
becomes: 


'' f 

a = — M I N dh. (5) 

f2 J 

The correct conversion of polarization rotation date to electron content date de- 
pends on the altitude chosen for the calculation of M. A method for arriving at 
such an altitude is shown in Figure 2. The height distribution of the electron 
density was obtained by converting the topside and bottomside ionograms recorded 
at close geographic proximity to Fort Monmouth into electron density profiles. 
Polarization- rotation date on 40 MHz was obtained from an overhead orbit of the 
S-66 satellite (nominal altitude = 100 km). The vertical component of the geo- 
magnetic field (curve B) was derived from the spikes in the topside ionogram 
which measure the gyrofrequency at the satellite altitude. At lower altitudes, 
the field intensity was calculated from the reference at th<^ satellite using the 
magnetic dipole approximation. Numerical, successive integration of the pro- 
duct N • B produced the curve of the polarization angle rotation E a which ends 
with the total polarization angle observed at the groimd ("measured*' in Figure 
2). The next procedure was integration of the bottomsMe and topside profiles 
to obtain / N dh. Next, the magnetic field component, M, was obtained from 
Equation 5 by inserting the measured total polarization angle Ea, and / N dh of 
the Integrated profiles. The figure indicates that for the above example M cor- 
responds to a local field value about 60km above hp^^^ . Figure 2 indicates fur- 
ther that 90% of the total polarization-angle rotation took place below 550 km. 


THE DATA 

The need for an empirical model of the global distribution of the ionospheric 
electron content for prediction purposes has focussed attention on the availability 
of such data on a global basis and on future requirements. The advent of beacon 
emissions from geostationary satellites has clearly facilitated date gathering at 
wide geographic areas, with the diurnal variation of content being obtained with- 
out contamination by satellite orbital changes. With low flying satellites, months 
of date from a particular station are necessary to obtain a diurnal variation. The 
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Figure 2. Calculation of ttie Faraday Polarization Rotatim Angle £a as a Func- 
tl<mofAltitudeUsingtheActuallyMea 80 redElectr(m Density Profile N(h), and the 
Geomagnetic Field Compcmcnt B Obtained from Gyrofrequency Spikes in the 
Tq;)side lonogram. 
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day-to-day as well as the geographic and seasonal variability of the ionosidiere 
is superimposed on the diurnal variation ttius oMained. 

Polarizati<» rc^tion measurements are performed continuously at Fort 
Monmouth, N. J. (40.25°N, 74.025°W) utilizing beacon emissions of the ATS-3 
(137.35 MHz). The subionospheric point (below 350km along the path from Fort 
Monmouth to ATS-3) is located at 37.1°N, 73.6*’W. For tbe purpose of this re- 
port, data observations for time periods betwe«i the 123rd and 178th days of 
1973, i.e. , May 3 to June 27, 1973, are presented. The data have been norma- 
lized to the vertical direction and, hence, represoit ttie vertical total electrcm 
cmitent. Representative dhimal variations are shown in Figures 3 throu^ 5. 

The following observations are made with respect to the absolute values of the 
ccmtent (and, equivalently, of the ionosi^eric time delay): 

1. rhe diurnal variati<m as well as the day-to-day variability of the total 
vertical content is evidrat. The lower and iqpper limits of the ccmtent 
are ~0.2 x 10*’ and ~3.0 x 10*’ respectively and correspondingly the 
i<mosi^ric time delay is -^1 nanosecond and 15.7 nanoseconds. Differ- 
ences of iq> to 100% can be observed in the figures. 

2. For file reported data, buildtqi of ionospheric ionization started daily at 
-'0400 EST. For file time of year considered, the dawn buildiqi begins 
when solar illumination is at —100 km (solar zraith angle = —98°). T<^ 
side i(»osi^eric densities have been reported to decrease with increasing 
solar illumination until about ground sunrise. ^ The discrepancy between 
topside density and total c<mtoit could be eiqplained fay thermally induced 
particle fluxes from the topside ionosphere to the bc^mside ionosphere 
during ionos|rfieric sunrise and by the quicker buildup of the bottomside 
i(Hiospheie. 

3. The buildup phase of the diurnal variation is smooth, and its time slope 
is nearly ccmstant. On most days, file sustained buildiq* phase ends prior 
to noontime, althou^ the diurnal maximum is not reached until later. 

4. The maximum of the content is variable in absolute value as well as in 
the time of occurrence. Such maxima usually occur after 1600 EST, and 
sometimes much later. 

5. In terms of absolute values of the electr<» content, magnetic activity does 
not have a ccmsistent effect on the content. Figures 3 fiirough 5 show 
cases when magnetically disturbed days show depressed ccmtent values 
with respect to quiet days, and vice versa. Even data taken on consecu- 
tive disturbed days have totally different characters. 
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Figure 3. Total Integrated Electron Content Diurnal Variation for the Days 123 Through 
129, 1973. (QQ, Q, D refer to Magnetically Extremely Quiet, Quiet, or Disturbed 
Days, Respectively). 


TOTAL ELECTRON CONTENT* K>' 


EXCESS TIME DELAYll.6GHi)_* 



TIME(EST) 


Figure 4. Same as 3, but for the Days 158 Through 163, 1973, 
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Figure 5. Same as 3, but for the Days 173 Through 178, 1973. 





6. Hie sustained decay of ionization with the setting of the sun starts, in 
most cases, from the maximum absolute value of the content. As during 
the buildup period, the decay with time is smooth and has a constant 
sl(^. Hie decay slope is sharper than the buildup slqpe. 

7. At nig^t the content continues to decline until it reaches an absolute mini- 
mim, from which the buildiqi phase starts at sunrise. The minimum may 
be reached just before sunrise or it may be sustained an hour or so before 
sunrise. 

8. Smaller scale variations are superimposed on the overall diurnal varia- 
tion of die electron content as is evident from Figures 3 through 5. These 
variati<ms are caused by ionospheric ionization irregularities along the 
signals' path. The irregularities are Hnite in extent and are Imown to be 
moving. Some of the irregularities are caused by ionospherically travel- 
ing internal gravity waves. ^ Their finite size m^ms that two relatively 
closely-spaced ground stations monitoring die same satellite may observe 
different vertically normalized electrcm content values since one's path 

to the satellite goes throu^ the irregularity while the other's does not. 

The traveling ionospheric disturbance uiiich cannot be predicted, except possibly 
statistically, constitutes a natural limit to the accuracy of any real-time predic- 
tion of the total integrated electron ccmtent. To ascertain the extent of this limit 
and its possible diurnal variation, the following procedure was used. Each day 
was divided into 6 four-hour intervals and the maximum deviation from the am- 
bient content at any of the intervals was recorded. The bounds of the deviation 
ranged from no deviation to an ~20% deviation. The deviations were then aver- 
aged for all the days and plotted in Figure 6. The following conclusions may be 
drawn: 

a. Electron content deviations occur, on the average, at all diurnal periods 
with a maximum of 5.4% and a minimum of 2. 8%. 

b. On the average, the deviations are smallest during the sustained iono- 
spheric buildup period (0400-0800 EST) and during the sustained iono- 
spheric decay period (1600-2000 EST). 

c. The greatest deviations occur during the day between the buildup and the 
decay. Nighttime deviations are also large. 

d. From the viewpoint of prediction schemes, the daytime deviations pose 
a greater problem to accuracy, since the daytime content is larger than 
the nighttime content by a factor of 3 or so. 
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Figure 6, Ionization Irregularity in Percent of Ambient Content 
for Six Daily Time Periods 


FUTITRE EXPERIMENTS 

The forthcoming launch of the geostationary ATS-F ^ will offer a unique and ex- 
cellent opportunity to study the diuinal behavior of the total electron content at 
various geographic locations. Our group will perform two experiments. In one 
experiment, the polarization rotation of the 140 MHz beacon will be measured. 
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This will not be different from the ei^riments carried out to date utilizing other 
geostationary satellites. In a second e^qieriment, continuous dispersive groi^) 
delay measurements will be made on 140 MHz and 360 MHz, both of which will 
be con^ilemented with ±1 MHz sideband emissions. With such a modulation on 
both frequencies, the distance between peaks of the modulation envelope is 300 
meters in free space and very little different in the ionosphere. The modulation 
mivel<^ travels at the group velocity and is retarded in phase by AP/300 cycles 
at the modidati<m frequmcy where 


AP = 



( 6 ) 


At ni^t (assuming /N ds 10*’ el/m ^), Ap ~200 meters for the 140 MHz fre- 

quency and AP 30 meters for the 360 MHz frequency, i.e. , phase reduction is 
<1 cycle for both frequencies. Hie modulation phase difference is thus obtained 
unambiguously and yields the absolute value of the total electron content. For 
other times of the day when the total content is larger, the absolute value could 
easily be derived from the ni^ttime values since the total content varies contin- 
uously with time. 


Since the phase effect is independent of the terrestrial magnetic field, it gives 
the true electron content along the entire propagation path, as opposed to the 
content given by the polarization rotation effect which pertains only to a signal 
path portion of 1500km from the earth's surface. Hie difference between simul- 
taneous polarization and groiqi-delay data will yield the currently unknown exo- 
spheric electron content, i.e., the content between 1500km and s 3 mchronous 
altitudes. 


SUMMARY 

With stringent precision requirements, a navigation system which measures 
signal propagation time between high-altitude satellites and users at or near the 
earth’s surface must take into account the excess time delay of the signal owing 
to the existence of free electrons along its path. The excess signal time delay 
is proportional to the total electron content along the propagation path. The up>- 
per bound for non-anomalous ionospheric conditions is ~150 nanoseconds for a 
highly oblique propagation path at a frequency near 1600 MHz. 

For dual-frequency user-capability, the excess time delay may be measured in 
real-time and a correction applied to the apparent range. For single-frequency 
user-capability, the ionospheric time delay must be arrived at through empirical 
global prediction techniques based on existing and future electron content data. 
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The pertinent information will be transmitted to the user via the navigating sig- 
nal. To date no truly global coverage of the total electron content exists. Most 
existing data are based c«i Faraday rotation measurements of the electron con- 
tent in the earth's inunediate vicinity <mly (i.e. , up to ^1500 km). Future ex- 
periments will yield content values up to geostati(mary altitudes. 

A limit to the accuracy of the prediction is set by superposition of the variations 
of ionization irregularities on the diurnal variation of the total electron content. 
The greatest average percent variation was observed during the day between the 
buildiq) and decay phases of the diumal variation and amounted to '^5.4% of the 
ambient total electron content. Since the ccmtent is at its highest during the day, 
such perr 3ntage deviaticms cause relatively hi^ absolute time delay deviations. 

Of course the complexity of the ionospheric processes and their interrelationships 
witii many otiier geophysical phenomena will introduce additional uncertainties 
into tte predictions. 
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QUESTION AND ANSWER PERIOD 


MR. EASTON: 

Any questions ? No questions ? 

(No response.) 

MR. EASTON: 

Thank you very much, Mr, Soicher. 
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